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ABSTRACT 

The  detection  performance  of  microwave  search  radars  operating  in  close  proximity  to 
the  sea  surface  is  evaluated.  The  effects  of  media  characteristics  on  the  propagation  of 
electromagnetic  waves  have  been  incorporated.  Specular  and  diffused  scattering  from  a 
rough  surface,  and  the  effects  of  the  curvature  of  the  earth's  surface  have  been  included 
in  the  study.  Additionally,  surface  ducting  effects  caused  by  atmospheric  anomalies  are 
presented.  Some  design  and  operational  considerations  which  can  improve  the  detection 
performance  of  a  surface  search  radar,  are  also  presented. 
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I.  INTRODUCTION 

A.  MOTIVATION  OF  STUDY 

Microwave  search  radars  installed  on  ships,  or  at  shore  stations,  operating  close 
to  the  sea  surface,  have  a  strong  interaction  with  the  propagating  medium  and  the  sea. 
The  reflected  field  is  a  function  of  the  transmitting  frequency,  the  incidence  angle,  and 
the  electric  properties  of  the  surface  boundary.  The  atmospheric  medium  also  interferes 
with  the  electromagnetic  wave  by  refraction,  absorption,  and  attenuation. 

The  present  study  tries  to  evaluate  quantitatively  the  effects  of  sea  roughness 
and  special  atmospheric  conditions,  like  ducting,  on  electromagnetic  propagation.  It  is 
shown  that  the  propagation  is  considerably  influenced  by  the  media  between  the  radar 
and  the  target. 

The  motivation  of  analyzing  the  irregularities  of  the  sea  surface  with  periodic  or 
random  variations,  together  with  atmospheric  refractivity,  is  to  predict  and  optimize 
the  detection  ranges  of  a  search  radar. 

B.  SCOPE  OF  WORK 

This  study  begins  by  introducing  the  Rayleigh  criterion  and  the  Fresnel  zone 
model  for  discriminating  between  a  rough  and  a  smooth  surface.  The  specularly  re- 
flected wave  interferes  with  the  direct  propagating  wave,  resulting  in  a  constructive 
and  destructive  interference. 

Computation  and  illustration  of  detection  patterns,  depending  on  radar  parame- 
ters and  sea  surface  characteristics,  are  performed.  The  resulting  plots  show  a  strong 
effect  of  constructive  and  destructive  interference  between  the  direct  and  reflected  fields. 


It  is  shown  that  under  some  special  conditions  the  detection  range  could  be  doubled 
relatively  to  free  space  range. 

The  abnormal  propagation  caused  by  strong  vertical  gradient  of  tropospheric 
refractivity,  results  in  a  considerably  extended  horizontal  detection  range.  From  con- 
servation of  energy  point  of  view,  the  increased  density  of  energy  through  one  path 
results  in  a  decreased  density  through  other  paths  causing  "holes"  in  detection. 

The  results  of  this  study  show  that  by  appropriate  design  and  operation,  detection 
ranges  of  a  radar  could  probably  be  improved.  Also,  the  limited  detection  range  could 
be  predicted  and  optimized. 

C.     ORDER  OF  WORK 

Chapter  II  presents  a  theoretical  study  on  surface  roughness,  specular  reflection 
and  diffused  scattering.  The  Rayleigh  criterion  and  Fresnel  zones  for  surface  roughness 
are  introduced. 

Chapter  III  introduces  a  simplified  model  for  propagation  above  a  spherical  sur- 
face. The  divergence  factor  and  the  reflection  coefficient  from  the  sea  surface  are 
discussed  in  detail.  The  shadowing  effect  for  radar  with  very  low  grazing  angles,  and 
the  depolarization  by  the  reflecting  surface  are  also  introduced. 

Chapter  IV  is  based  on  the  study  of  Chapter  II  and  III.  This  chapter  discusses 
the  lobing  phenomena.  Computer  plots  of  detection  contours  are  illustrated.  Conclu- 
sions regarding  improved  detection  range  for  low  flying  aircraft  and  surface  targets  are 
suggested. 

Chapter  V  introduces  abnormal  propagation  in  the  atmosphere.  The  atmospheric 
conditions  for  this  phenomena  to  occur,  and  the  effect  on  radar  detection  ranges  are 
discussed. 


Chapter  VI  is  the  summary  and  conclusion  of  this  study.  The  Appendix  gives  an 
analytical  derivation  of  Fresnel  zones. 


II.  RAYLEIGH  CRITERION  AND  FRESNEL 
ZONES  FOR  SURFACE  ROUGHNESS 

A.  INTRODUCTION 

When  an  electromagnetic  short  wave  (VHF  and  above)  propagates  above  the 
earth's  surface,  it  is  affected  by  the  earth's  surface  and  its  atmosphere.  The  reflected 
field  of  the  electromagnetic  wave  from  such  boundaries  depends  on  the  angle  of  in- 
cidence, the  wavelength  and  the  characteristics  of  the  reflecting  plane  (conductivity, 
permittivity,  and  permeability).  The  difficulty  of  qualitatively  analyzing  the  influence 
of  the  medium  on  the  propagating  wave  is  due  to  the  fluctuating  irregularities  of  the 
surface  boundaries.  These  irregularities  impose  limitations  on  the  performance  of  com- 
munication and  radar  systems.  These  random  fluctuations  of  the  parameters,  defining 
a  practical  terrain  profile,  are  the  main  obstacle  in  establishing  its  relationship  with 
the  properties  of  the  electromagnetic  wave.  For  this  reason  the  prediction  of  radar 
and  communication  performance  is  limited.  In  other  words,  the  effects  on  propagation 
introduces  degradations  in  radar  performance. 

In  this  chapter  we  will  discuss  the  theory  of  specular  reflection  and  its  effect 
upon  field  strength  on  the  assumption  that  the  reflection  coefficient  of  a  perfectly 
smooth  earth  are  known.  Effects  of  surface  roughness  will  also  be  discussed. 

B.  ANALYTICAL  MODEL  AND  CONSIDERATIONS 
1.       Surface  Roughness:  General  Definition 

The  first  stage  in  constructing  a  model  for  the  propagation  problem  is  to 
analyze  and  understand  the  reflecting  boundary  which  is  usually  varying  in  position 
and  time,  and  therefore  is  not  plane. 


The  basic  difference  between  a  rough  and  a  smooth  surface  is  that  a  smooth 
surface  will  reflect  the  incident  plane  wave  specularly  in  a  single  direction,  whilst  a 
rough  surface  will  scatter  the  wave  in  various  directions  (see  Fig  2.1,  2.2).  From  this 
definition  one  can  conclude  that  the  roughness  of  the  surface  is  determined  basically 
by  the  angle  of  incidence  and  the  frequency  of  the  plane  wave.  A  smooth  surface  is  a 
singular  limiting  case  of  a  rough  one. 


•  '///  f//s 
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Figure  2.1:  Reflected  Waves  from  Perfectly  Smooth  Surface  [Ref.  1]. 


In  Chapters  III  and  IV  we  will  introduce  the  sea  medium  which  will  include 
the  time  varying  nature  of  the  scattering  surface  and  its  effects  such  as  the  fading  of 
the  received  signal  by  the  shape  and  movement  of  the  water  waves  [Ref.  3]. 
2.     Fundamental  assumptions 

The  empirical  and  analytical  theory  of  rough  surface  scatter  is  usually  sim- 
plified by  some  of  the  following  assumptions  [after  Ref.  3]: 
1)   Far  field  approximations  only. 


Figure  2.2:  Reflected  Waves  from  Rough  Surface  [Ref.  2]. 


2)  Shadowing  effects  are  neglected. 

3)  The  radius  of  curvature  of  the  scattering  surface  is  much  larger  than  the  wave- 
length of  the  incident  wave. 

4)  Multiple  scattering  is  neglected. 

5)  Particular  models  of  surface  roughness  are  restricted  (periodic  and  random). 
3.     Roughness  Criteria 

As  mentioned  earlier,  the  efficiency  of  specular  reflection  from  a  rough 
surface  is  a  function  of  the  ratio  of  the  surface  roughness  dimensions  to  the  wavelength, 
and  of  the  grazing  angle.  We  will  try  to  define  the  values  of  the  grazing  angle,  the 
wavelength,  and  the  surface  roughness  at  which  the  specular  reflection  changes  into 
diffuse  scattering. 

The  best  known  criterion  involving  these  parameters  is  that  of  Lord  Rayleigh. 
The  path  difference  between  the  two  waves  is 

6  =  2AhsinB.  (2.1) 

The  phase  difference  A<f>  corresponding  to  this  path  difference  is: 

AwAh 
A<t>  =  k6  =  — —  sin0,  (2.2) 

A 

where  A  is  the  wavelength  and  the  propagation  number  is  given  by  k  =  y. 

For  small  phase  differences,  there  is  almost  no  effect  by  the  roughness  and 
the  surface  behaves  as  a  smooth  surface.  This  can  be  achieved  if  —  or  9  are  very  small 
which  leads  to  A<f>  «  0.  Increasing  —  and/or  0  will  increase  the  phase  difference. 
When  the  phase  difference  is  equal  to  n  the  two  waves  will  be  in  phase  opposition 
and  cancel  in  that  direction.  According  to  conservation  of  energy,  the  energy  of  the 
electromagnetic  wave  is  scattered  in  other  directions.  Thus  we  conclude  that  for: 

1)  A<j)  =  0  the  surface  reflects  specularly  and  is  smooth. 

2)  A(f)  =  7r  the  surface  scatters  and  is  rough  [Ref.  1]. 

A  criterion  for  discriminating  between  rough  and  smooth  surfaces  is  chosen 
arbitrarily  as  7r/2  which  is  the  average  phase  difference  between  these  two  extremes 


[Ref.  5].  Thus  a  surface  is  considered  to  be  smooth  for  the  height  of  roughness  Ah  to 
be  smaller  than: 

osin  a 
This  criterion,  sometimes  called  the  Rayleigh  criterion,  is  a  simplified  argument  based 
on  ray  theory  for  describing  the  irregularities  of  a  surface. 

The  Rayleigh  model  excludes,  by  definition  (because  of  the  ray  assumption), 
the  case  when  the  roughness  of  the  surface  elements  are  comparable  in  size  to  the 
wavelength  of  the  propagating  wave. 

The  sea  surface  in  such  a  case  usually  consists  of  irregular  distributions  of 
sizes  and  shapes  of  the  waves.  An  accurate  solution  in  this  case  requires  complex 
methods  not  yet  developed  to  satisfy  boundary  conditions  of  a  random  rough  surface. 
However,  in  the  literature,  simplifications  have  been  made  in  order  to  calculate  the 
effective  reflection  coefficient,  using  sea  surface  models  with  sinusoidal  or  wedge  shaped 
waves  [Ref.  3]. 

4.      Active  Scattering  Region 

As  mentioned  above,  the  complexity  of  the  real  surface  is  so  severe  that  no 
analytical  model  can  describe  it  accurately.  In  this  section  we  will  try  to  consider  the 
effects  of  scattering  in  terms  of  a  simple  model  called  the  Fresnel  zones,  or  the  glistening 
surface.  Ray  analysis  suggests  that  reflections  would  be  confined  to  a  small  region  of 
the  rough  surface,  which  depends  on  the  transmitter  receiver  configuration. 

The  glistening  surface,  or  the  Fresnel  near  zones,  is  defined  as  that  part  of 
a  rough  surface  that  reflects  a  significant  amount  of  electric  field  of  the  propagating 
wave  into  the  receiver,  at  given  heights  of  the  source  and  the  receiver  [Ref.  4]. 

Applying  Huygen's  principle  to  the  ray  theory,  one  should  interpret  the  re- 
flected ray  from  a  single  point  on  the  surface  as  the  fact  that  an  energy  source  illumi- 
nating a  large  portion  of  the  surface  creates  currents  which  radiate  in  all  directions.  At 
any  point  in  space  the  total  field  is  the  sum  of  the  radiation  from  these  currents  and 
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Figure  2.3:  Glistening  Zones  on  a  Reflecting  Plane  [Ref.  4]. 

the  direct  propagating  field.  Figure  2.3  illustrates  a  scenario  consisting  of  a  transmitter 
at  height  Z\  and  a  receiver  at  height  z2  operating  above  a  smooth  plane. 

The  direct  ray  path  is  Rd  and  the  indirect  path  from  a  point  on  earth  is 
Ri  +  R2.  As  mentioned,  the  incident  field  induces  currents  on  the  plane  xy  which  radiate 
a  reflected  field.  We  want  to  find  the  zones  in  which  the  path  difference  between  the 
direct  and  the  indirect  path  is  constant  and  equals  to  A<f>.  The  path  difference  is  [Ref. 
1,3]: 

Rl  +  R2-Rd  =  6.  (2.4) 

Since  the  direct  path  is  unchanged  and  S  is  constant: 

R1  +  R2  =  Const.  =  Rd  +  6  (2.5) 

Putting  the  path  difference  S  into  phase  difference  notation: 

A<f>  =  6-k,  (2.6) 


or 


A<f>  = 


6-2n 


(2.7) 


Equation  2.5  is  an  equation  of  ellipsoids  with  T  and  R  as  foci  and  the  direct  path 
Rd  as  axis  of  revolution  (see  Figure  2.3).  The  geometric  locus  of  all  points  on  the  xy 
plane  which  intersect  with  the  ellipsoids  are  ellipses  which  describe  a  constant  phase 
difference,  A0.  In  order  to  describe  more  zones  of  other  constant  path  difference  we 
increase  8,  for  instance,  by  A/2,  or  tt  radians  (for  phase  interpretations). 

We  find  a  family  of  ellipses  which  describes  zones  of  constant  phase.  Each 
zone  differs  from  its  adjacent  zone  by  7r,  so  the  radiation  from  two  successive  zones 
is  in  opposite  phase.  They  do  not  cancel  totally  because  of  the  slow  variation  in  the 
amplitude  between  zones.  We  shall  see  that  the  first  zone,  n  =  1,  is  the  contributing 
zone  [Ref.  4]. 

a.      Fresnel  zones  dimensions 

In  order  to  apply  smooth  earth  phase  difference  formulas  one  should 
know  the  dimensions  of  the  Fresnel  zones,  which  can  give  an  indication  as  to  what  the 
minimum  dimensions  of  a  smooth  area  should  be  for  these  formulas  to  be  valid. 

The  approximate  path  difference  between  the  direct  path  and  the  geo- 
metrical point  (only  one)  of  reflection  is: 

2z\Z2 


6o 


(2.8) 


the  nth  ellipse  path  difference  is  determined  by: 

6n&60  +  — . 


(2.9) 


By  using  this  approximation  and  substituting  it  into  Equations  A.  13, 

A.  15,  and  A.  11  and  using  the  assumptions  of  Equation  A.  16  (Appendix  A),  we  find 

useful  formulae  for  the  semiminor  axis  (yin),  the  semimajor  axis  (xi„)  and  the  center 

(x0n)  of  the  nth  ellipse: 

yn\r 


%\n 


Vln  «  ±- 


r 

X0n  ±   ~ 


2     ^1  +  ^ 


1  + 


nXr 


xln  —  x0n  i  £m 


(2.10) 
(2-11) 
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and 

*°»  * 2 ; ,  sfe  •  <2-12) 

'  nAr 

The  axea  of  the  nth  ellipse  is: 

An  =  7r£nyln.  (2.13) 

By  substituting  Equations  2.10  and  2.11  into  Equation  2.13,  we  obtain: 

4(l  +  ^£)^"  (        } 

b.     Special  Practical  Cases 

The  minimum  area  for  smooth  surface  assumption  (Fresnel  zone)  calcu- 
lated by  Equation  2.14  can  be  simplified  for  special  practical  cases  of  a  radar  tracking 
an  aircraft: 

^  «  1.  (2.15) 

22 

The  height  of  the  radar  antenna  is  z\  and  the  height  of  an  aircraft  is 
z2.    Since  the  first  zone  is  the  most  contributing  area  to  the  total  field  received  at  R 
(Fig.  2.3),  we  should  consider  only  this  zone  (n  =  1).   At  the  reflection  point  for  the 
lowest  lobe  of  the  interference  pattern,  r  ^>  Zi,z2  and  S  =  ^. 
The  parameters  yin,#in  and  z0n  are  simplified  to  [Ref.  4]: 

y„    a    2v/2*i«T^r,  (2.16) 


2\Z2z-i 

Xq 

3n/2zi   ' 
A 

Xoi  ± 

A     ' 

122? 

xu    «    x0i±— y-^,  (2.17) 

*oi    «    -y1.  (2.18) 

Two  cases  are  considered  to  illustrate  the  dimensions  and  positions  of 
the  first  Fresnel  zone.  In  the  first  case,  the  range,  the  height  of  the  source  z\ ,  and  the 
height  of  the  receiver  z2,  are  given  and  it  is  required  to  determine  the  position  and 
dimensions  of  the  first  zone.  In  the  second  case  r  and  z\  are  given  and  we  want  to  find 
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dimensions  of  the  first  zone.  In  the  second  case  r  and  Z\  are  given  and  we  want  to  find 
the  dimensions  and  position  of  the  first  zone  which  will  produce  the  reflected  wave  to 
form  the  lowest  lobe  of  the  interference  pattern  for  which  S0  =  £.  The  calculated  values 
for  given  grazing  angle  and  mean  height  of  sea  waves  are  given  in  Table  2.1  [Ref.  3,4]. 

TABLE  2.1:  FRESNEL  ZONES  AREA  Ax  [REF.  4]. 


Parameter 

Case  1 

Case  2 

A  [cm] 

10 

100 

10 

100 

r  [m] 

10,000 

10,000 

10,000 

10,000 

z2  [m] 

1000 

1000 

2.5 

2.5 

Z\  [m] 

100 

100 

100 

100 

4, 

20m 

20m 

5cm 

50cm 

x0i  [m] 

910 

943 

9350 

6830 

yn  [m] 

9.1 

29 

6.6 

44 

Mm] 

83 

262 

615 

2760 

*[°] 

6.3° 

6.3° 

0.59° 

0.72° 

frmax   [m] 

0.11 

1.1 

1.22 

1.25 

A,  [m2] 

2,370 

23,900 

128,000 

382,000 

Table  2.1  indicates  that  the  minimum  areas  of  the  first  zone  for  smooth 
assumptions  seem  to  be  large,  and  are  larger  with  increasing  wavelength.  For  mean 
peak  to  peak  wavelength  of  almost  lm,  in  case  2,  the  smooth  area  for  the  first  lobe  of 
interference  is  much  larger  than  the  smooth  area  in  case  1,  especially  for  short  wave 
lengths. 

The  contribution  of  energy  to  the  receiver  from  each  zone  is  calculated 
by  integrating  the  current  excitation  of  each  zone.  Though  the  successive  zones  differ 
in  phase  by  n  radians,  it  is  not  evident  that  their  contribution  to  the  received  energy 
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cancels  out  totally  and  that  only  the  first  zone  or  part  of  it  contributes  to  the  received 
energy. 

The  reason  for  this  doubt  could  be  caused  by  the  variation  in  the  am- 
plitude of  excitation  within  a  single  zone,  especially  in  small  grazing  angles.  From  later 
research  it  seems  that  the  first  Fresnel  zone  is  likely  to  be  the  most  important  zone  of 
contribution  to  the  received  energy,  together  with  two  other  areas  near  the  source  and 
near  the  receiver  which  differ  from  the  Fresnel  zones. 

In  this  discussion  only  the  phase  of  the  reflected  radiation  has  been  con- 
sidered. We  did  not  consider  the  conductivity  and  dielectric  constant  nor  the  reflection 
coefficient.  These  will  be  discussed  in  the  next  chapter.  The  sea  surface  is  more  com- 
plex than  the  land  surface,  with  time  variations  in  addition  to  fluctuations  in  amplitude 
and  phase  of  the  scatterers. 
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III.  EARTH'S  CHARACTERISTICS 
AFFECTING  PROPAGATION 

A.  INTRODUCTION 

After  discussing  the  roughness  of  the  surface  in  Chapter  II  and  deriving  roughness 
criteria  and  a  model,  we  want  to  discuss  in  this  chapter  the  major  characteristics  of  the 
earth's  surface  which  will  affect  the  magnitude,  phase,  and  direction  of  a  propagating 
wave. 

Firstly,  we  will  consider  the  effects  of  a  spherical  surface  and  refraction  effects  and 
present  a  model  for  horizon  propagation  range,  taking  into  account  the  convex  shape 
of  the  earth's  surface.  Secondly  we  will  analyze  the  reflectivity  of  the  surface. 

Compared  with  propagation  in  free  space,  the  presence  of  terrain  affects  the  Radar 
Cross  Section  (RCS)  of  scatterers.  The  reflection  from  terrain  creates  an  interference 
pattern  which  will  be  discussed  in  detail  in  Chapter  IV.  Also,  the  terrain  casts  a  shadow 
which  introduces  diffraction  phenomena.  The  main  outcome  will  be  a  reflection  coef- 
ficient for  a  rough  and  a  smooth  surface,  which  will  define  quantitatively  the  reflected 
and/or  the  scattered  electromagnetic  wave;  especially  from  the  sea-surface. 

Finally  the  shadowing  and  depolarization  effects  will  be  introduced.  These  effects 
are  also  significant  in  determining  the  reflected  and/or  scattered  field  of  an  electromag- 
netic propagating  wave. 

B.  SPHERICAL  SURFACE  AND  REFRACTION  EFFECTS 

Electromagnetic  wave  reflectivity  becomes  more  complex  when  the  surface  is  not 
flat.   The  reflected  wave  is  decreased  in  amplitude;  also  the  phase  difference  between 
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propagation  over  a  convex  surface.  The  difference  between  the  plane  and  spherical 
surface  is  that  in  spherical  surface  the  height  of  the  transmit  and  receive  antennas  is 
increased  in  order  to  achieve  equivalent  line  of  sight.  Figure  3.2  illustrates  the  reduction 
in  the  reflected  wave's  amplitude  towards  a  specific  direction  due  to  the  divergence  of 
the  waves  into  a  larger  area,  compared  with  a  flat  surface.  [Ref.  2] 

Realizing  that  the  propagation  of  the  electromagnetic  wave  is  in  the  atmosphere, 
one  must  consider  the  refraction  effect.  This  effect  causes  the  waves  to  bend  during 
propagation  rather  than  propagating  in  straight  lines,  usually  resulting  in  an  extension 
of  the  maximum  range  of  propagation.  It  is  convenient  to  analyze  the  effect  of  the 
phenomena  together  with  the  earth's  curvature.  Detailed  discussion  on  atmospheric 
refraction  is  carried  out  in  Chapter  V. 

1.     Horizon  Range 

The  reduced  heights  (h[,  h'2)  of  the  radar  and  target  for  a  plain  surface  are 
calculated  by: 

*-*,-£,  (3.1) 


and 


cP 

ho  =  hi — ,  (3.2) 

2*re' 


where  dx  is  the  curved  radar  horizon  when  h[  is  zero,  d2  is  the  curved  distance  from 
the  reflecting  point  to  the  receiving  antenna  and  hi  and  h2  are  the  actual  heights  of 
the  transmitting  and  receiving  antennas  respectively  (see  Figure  3.1).  The  average 
effect  of  standard  atmospheric  refraction  is  interpreted  by  increasing  the  earth  radius 
rc  by  the  factor  k  =  4/3.  (For  more  details  see  Chapter  V.)  The  differences  in  heights 
of  the  radar  antenna  and  the  target,  due  to  refraction  and  earth  curvature,  are  given 
by  [Ref.  2]: 
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Figure  3.1:  Spherical  Surface  Geometry  [Ref.  2], 
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Figure  3.2:  Divergence  by  a  Convex  Surface  [Ref.  2], 
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Afclt2  =  -f.  (3.3) 

The  radar  horizon  dh  is  given  by: 

dh  =  ^2krehi  +  h\  +  y/2kreh2  +  /*2-  (3-4) 

The  distance  c^  +  d2  is  the  horizon  range  of  the  radar  if  h2  is  zero.    Using 
this  fact  in  Equation  3.3,  the  horizon  range  then  can  be  approximated  to: 


ih  w 


\J2kreh,  (3.5) 


where  h  is  the  antenna  height  of  the  radar  and  kre  is  the  effective  earth  radius.  For 
standard  atmospheric  conditions  (k  =  4/3)  the  distance  dh  has  the  form  of: 

dh  =  1.23v//i(ft)  in  nautical  miles.  (3.6) 

Furthermore,  the  relationship  between  the  optical,  radar  and  geometric  (prop- 
agating along  perfect  straight  lines)  horizons  for  standard  atmosphere  is:  [Ref.  8] 

dh  =  1.074  =  1.15<f5,  (3.7) 

where  dh,  d0  and  dg  are  the  radar,  optical,  and  geometrical  horizon  ranges.  Equation  3.7 
stresses  that  under  standard  conditions  the  radar  horizon  range  is  7%  larger  than  the 
optical  horizon  and  15%  larger  than  the  geometric  horizon.  In  abnormal  conditions  such 
as  ducting,  the  electromagnetic  waves  can  propagate  to  longer  ranges  due  to  diffraction 
from  the  earth,  combined  with  refraction  from  atmospheric  layers.  [Ref.  2] 

C.      COMBINED  REFLECTION  COEFFICIENT 

As  mentioned  above  the  physical  principle  which  allows  radar  coverage  to  be 
extended  beyond  the  geometrical  horizon  is  divergence  or  refraction.  The  electromag- 
netic wave  is  diffracted  around  the  earth's  curvature,  the  range  of  propagation  depends 
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on  the  frequency.  The  larger  the  wavelength  (lower  frequency)  the  more  the  wave  is 
propagated,  due  to  refraction. 

This  phenomena  is  applied  in  VLF  (Very  Low  Frequency)  for  world-wide  com- 
munications. Since  radar  frequencies  have  short  wavelengths  compared  with  earth 
dimensions,  only  small  amounts  of  energy  will  be  diffracted,  resulting  in  insufficient 
propagation  beyond  the  radar  horizon  by  divergence.  [Ref.  5] 

The  combined  reflection  coefficient  of  earth  is  determined  by  three  factors: 

1)  D,  the  divergence  factor  that  describes  the  reduction  in  reflection  due  to  earth's 
curvature  (see  Figure  3.2), 

2)  p,  the  reflection  coefficient  of  a  plane  smooth  surface,  and, 

3)  R,  a  factor  which  defines  the  surface  roughness. 

The  combined  reflection  coefficient  is  the  product  of  the  three  factors  mentioned 
above  [Ref.  2]: 

combined 

reflection     »  =  pDR.  (3.8) 

coefficient 

The  factor  p  may  be  calculated  from  Fresnel's  equations  (see  next  section).  The 
factor  R  is  equal  to  unity  if  the  surface  is  smooth.  If  the  surface  is  extremely  rough  R 
is  equal  to  zero.  The  factor  D  varies  between  0  and  1,  and  D  approaches  unity  when: 

1)  hi  or  h2  (Figure  3.1)  approaches  zero  or  k  becomes  very  large  (especially  over 
water),  and, 

2)  For  radar  geometry  such  that  the  earth  can  be  assumed  flat  [Ref.  2], 
The  divergence  coefficient  D  is  given  by: 

-1/2 


D  ~ 


1  + 


2did,2 


kredh  sin  6 


(3.9) 


where  9  is  the  grazing  angle. 
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D.     REFLECTION  COEFFICIENT  OF  PLANE  SMOOTH  EARTH 

Prior  to  discussing  the  characteristics  of  reflection  of  a  plane  wave  from  rough 
surface,  one  should  first  consider  the  reflection  characteristics  from  a  smooth  surface. 
The  relative  complex  dielectric  constant  erc  is  given  by: 

erc  =  —  -  6(h'A<7,  (3.10) 

where  e0  is  the  dielectric  constant  of  free  space,  A  is  the  wavelength  in  meters,  a  is  the 
conductivity  in  mho  —  m-1)  and  e  is  the  dielectric  constant. 
The  relative  magnetic  permeability  is  given  by  [Ref.  1]. 

Hr  =  — ,  (3.H) 

where  [i  is  the  magnetic  permeability  of  earth  and  fx0  is  that  of  free  space.  As  mentioned 
in  Section  C,  the  reflection  coefficient  of  a  plane  smooth  surface  can  be  calculated  by 
Fresnel's  equations  [Ref.  2,3]: 

_  sina' —  (erc  —  cos2a')1/2  ,,  _  „_ 

r"  =  ~ TT^ r^UT  =  PH<r*H,  (3-12) 

sin  a  +  (erc  —  cos'  a)1' l 

r         ercsma'  -  (erc  -  cos2  a')1'2  _j4> 

erc  sin  a  +  (erc  —  cos'  a  )1/z 

where  a'  denotes  the  incident  and  reflected  angle  from  a  tilted  surface  (see  Figure  3.3), 
and  4>h  and  (f>r  are  the  horizontal  and  vertical  phase  of  the  reflected  wave.  The  mag- 
nitude and  phase  of  the  reflection  coefficient  depends  on  the  polarization,  frequency 
and  angle  of  incidence  of  the  electromagnetic  wave;  and  on  the  dielectric  constant  and 
conductivity  of  the  surface. 

The  Fresnel  Equations  3.12  and  3.13  correspond  to  two  polarizations  of  the  electric 
field.  The  subscripts  H  and  V  stand  for  horizontal  and  vertical  polarizations  respec- 
tively. Fresnel  equations  are  given  under  the  assumption  that  the  wave  is  propagating 
in  free  space,  and  that  the  dielectric  surface  is  non-magnetic. 
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When  a  is  equal  to  |  (normal  incidence),  the  horizontal  and  vertical  reflection 
coefficients  are  identical  in  magnitude  but  opposite  in  sign.  In  this  case  Equations  3.12 
and  3.13  become: 

1  -  €X'2 

r*  =  —-172,  (3-14) 

1  +€rc 

e    -  e1/2 
Tv  =  ":£—^  =  -IV  (3.15) 

E.     ILLUSTRATIONS  OF  REFLECTION  COEFFICIENT  OF  A 
SMOOTH  SURFACE 

Figures  3.4  to  3.7  illustrate  the  dependence  of  the  magnitude,  p,  and  the  phase, 
(f>;  on  the  grazing  angle  9,  wavelength  A,  and  the  polarization.  In  the  case  of  vertical 
polarization  there  is  a  sharp  decrease  in  the  reflection  coefficient  due  to  the  increase 
of  transmission  into  the  surface  when  the  angle  of  incidence  is  close  to  Brewster  angle, 
(as).  At  this  angle  py  reaches  its  minimum.  The  phase  (f>v  is  [Ref.  2,3]: 

<t>v  =  -  -aB.  (3.16) 

In  Figure  3.6  it  is  shown  that  <f>v  depends  also  on  the  wavelength  of  electromagnetic 
wave. 

When  the  electric  field  is  horizontally  polarized,  only  slight  variations  occur  in  the 
magnitude  and  phase  of  the  reflection  coefficient  (pn  and  (f>n  respectively),  with  respect 
to  the  incidence  angle  (See  Figures  3.4  and  3.5).  The  phase  <f>jj  is  almost  180°  for  long 
wavelengths  and  increases  slightly  with  incidence  angle.  In  many  practical  applications 
Ph  and  <f>H  are  approximately  unity  and  180  degrees  (7r  radians)  respectively,  i.e., 
T  =  -1.  [Ref.  1,2] 

For  very  small  incident  angles  py  and  <f>v  are  also  close  to  unity  and  7r,  respec- 
tively. Such  a  case  occurs  for  anti-ship  sea-skimming  missiles  guided  by  a  radar  and  for 
microwave  searching  radars  operating  close  to  the  ocean  surface.  For  incidence  angles 
larger  than  10°,  <f>y  decays  to  zero  (see  Figure  3.5).    The  characteristics  mentioned 
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Figure  3.3:    Incident  and  Reflected  Wave  by  a  Smooth  Dielectric  Surface 
[From  Ref.  2]. 
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above  for  small  incidence  or  grazing  angles  (<  1°)  will  play  a  major  role  also  when  the 
surface  of  incidence  is  rough.  A  complex  effect  of  shadowing  arises  for  small  incidence 
angles  and  rough  surfaces.  This  effect  will  be  discussed  in  the  next  section. 
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Figure  3.4:  Reflection  Coefficient  of  Very  Smooth  Sea  Water  [Ref.  2]. 

F.     THE  EFFECT  OF  SURFACE  ROUGHNESS  ON  THE  REFLECTION 
OF  ELECTROMAGNETIC  WAVES 

In  Chapter  II  we  discussed  the  Rayleigh  criterion  and  Fresnel  zones  for  smooth 
surface  approximation.  In  this  section  we  will  try  to  consider  some  more  aspects  of 
distinguishing  between  a  rough  and  a  smooth  surface.  Using  Equation  2.4  for  Rayleigh's 
criterion  of  smooth  surface  we  get  [Ref.  2,3]: 


AhsinO       1 

-a—  Ki- 


(3.17) 


From  this  simplified  criterion  one  would  expect  the  reflection  coefficient  to  be  a 
function  of  wavelength,  surface  roughness,  and  grazing  angle.  This  criterion  gives  only 
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Figure  3.5:  Phase  of  Reflection  Coefficient  <f>  of  Very  Smooth  Sea  [Ref.  2]. 


Figure  3.6:  Reflection  Coefficient  for  Low  Depression  Angles  of  Very  Smooth 
Sea  [Ref.  2]. 
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Figure  3.7:    Phase  of  Reflection  Coefficient  for  Low  Depression  Angle  of 
Very  Smooth  Sea  [Ref.  2], 

an  approximate  indication,  dealing  only  with  orders  of  magnitude.  From  empirical  data 
it  had  been  indicated  that  for  large  values  of  the  parameter  in  Equation  3.17,  where: 


Ah  sin  0        1 


A 


>V 


(3.18) 


A/isin0 


the  reflection  coefficient  is  independent  of  the  parameter  — jp-  for  a  specific  surface. 
However,  small  values  of  — ^-^  the  reflection  coefficient  increases  with  an  increase 
of  the  parameter.  Actually  there  is  no  sharp  distinction  between  rough  and  smooth 
surface.  [Ref.  4] 

One  fundamental  result  from  analyzing  the  scattering  field  from  a  rough  surface  is 
that  it  consists  of  two  components,  a  coherent  or  specular  component  and  an  incoherent 
or  diffused  component.  [Ref.  1] 

The  specular  reflection  is  defined  as  a  single  reflection  of  electromagnetic  wave 
from  a  rough  surface  in  the  direction  of  the  receiver.  Its  reflection  obeys  the  optical 
laws  and  its  phase  is  coherent.  In  other  words,  the  reflected  waves  from  a  rough  surface 
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can  be  considered  as  specular  if  they  arrive  at  the  receiver  with  equal  phases  and  have 
relatively  small  amplitude  fluctuations.  These  reflections  axe  from  the  Fresnel  zones  as 
described  in  Chapter  II. 

The  diffused  scattering  is  scattered  randomly  from  a  much  larger  area  than  the 
first  Fresnel  zone.  Its  phase  is  uniformly  distributed,  but  the  large  fluctuations  of  its 
amplitude  are  Rayleigh  distributed.  [Ref.  1] 

These  two  components  were  observed  in  several  experiments.  The  specular  re- 
flected field,  when  combined  with  the  direct  field,  will  produce  at  a  test  receiver  a 
periodically  varying  field.  When  a  diffuse  field  is  added  to  the  direct  field,  the  result 
will  be  a  randomly  varying  field.  [Ref.  1,3] 

Specular  reflection  and  diffuse  scattering  often  occur  simultaneously.  The  diffused 
scattered  field  could  be  the  reflection  from  large  numbers  of  uncorrelated  scatterers. 
It  would  be  expected  that  the  reflection  is  totally  diffused  and  independent  of  phase 
difference  A0.  This  is  not  always  true,  since  in  the  case  of  perfectly  smooth  surfaces 
the  diffused  field  is  zero  [Ref. 4].  Also  the  specular  field  many  be  produced  by  many 
scatterers  since  it  can  be  produced  by  a  linear  summation  of  many  specular  fields.  In 
reality,  the  total  field  is  the  sum  of  both  the  specular  and  diffused  fields. 

Many  attempts  were  made  to  model  a  reliable  specular  reflection  coefficient, 
caused  by  a  Gaussian  distributed  surface  as  a  function  of  mean  wave  height  of  the 
sea,  grazing  or  incident  angle  6,  and  the  transmitting  wavelength  A.  The  improved  re- 
flection coefficient  pa  does  not  take  into  account  the  shadowing  effect.  The  magnitude 
of  the  specular  reflection  coefficient  is  calculated  by  [Ref.  2]: 

p.  =  po  ■  DR.,  (3.19) 

where  R3  is  the  specular  scattering  coefficient  for  rough  surface,  po  is  the  reflection 
coefficient  of  a  plane  and  perfectly  smooth  surface,  and  D  is  the  divergence  coefficient 
due  to  earth's  curvature.  Modifications  in  the  reflection  coefficient  caused  by  a  Gaussian 
distributed  rough  surface  was  calculated  by  Beckman  [Ref.     1].     The  mean  square 
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variations  in  Rs  is  given  by: 

<R2>=e~2W.  (3.20) 

The  phase  difference  A<f>  is  calculated  by: 

A<j>  =  27rAhsm0/\,  (3.21) 

where  Ah  represents  the  standard  deviation  of  the  Gaussian  distribution  of  wave  heights 
[Ref.  2].  Equation  3.21  is  an  analogy  of  phase  difference  for  Rayleigh  criterion  (Equa- 
tions 2.2,  2.4).  The  improvement  in  Equation  3.21  is  that  Ah  is  not  just  a  magnitude 
but  defined  more  correctly. 

Most  of  the  properties  of  reflection  from  a  smooth  surface  can  be  applied  to 
specular  reflection  from  a  rough  surface.  The  basic  difference  is  that  specular  reflection 
[Ref.  1]: 

1)  fluctuates,  and, 

2)  is  smaller  in  amplitude. 

These  facts  are  taken  into  account  through  the  specular  scattering  coefficient  Rs. 

Different  models  for  <  R2  >  were  suggested.  In  all  of  these  <  R2a  >  is  a  function 
of  the  phase  difference  A<f).  The  additional  models  of  <  R2  >  are  [Ref.  1]: 

1)  a  uniformly  distributed  surface: 

<  R]  >=  sinc2A<£,  (3.22) 

2)  first  Fresnel  zone  of  a  uniformly  distributed  surface: 

<  R]  >=  (1  -  A<f>/Tr)smc2A(f>.  (3.23) 

By  checking  these  three  models  (Equations  3.20,  3.22,  3.23)  against  experimental 
results,  the  improved  Gaussian  surface  distribution,  Equation  3.20,  gave  the  closest 
results  to  the  experimental  ones  [Ref.  2].  Figure  3.8  illustrates  the  dependence  of 
the  mean  square  of  the  reflection  coefficient,  on  the  phase  variance,  for  the  improved 
Gaussian  model. 
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From  experiments  it  was  seen  that  there  were  some  variations  of  <  R?a  >  from  the 
curve  when  the  abscissa  is  >  120  mradians.  Those  variations  were  caused,  probably,  by 
diffused  scattering;  or  by  the  combination  of  diffused  and  specular,  scattering/reflection 
[Ref.  2]. 

Until  now  we  have  discussed,  quantitatively,  only  the  specular  reflection  from 
a  rough  surface.  Now  we  want  to  define  and  analyze  the  diffused  scattering.  The 
reflection  coefficient  of  diffuse  scattering  pd  can  be  written  similarly  to  the  specular 
reflection  coefficient: 

pd  =  poDRj,  (3.24) 

where  Rj.  is  the  coefficient  which  defines  the  irregularities  of  the  surface,  and  po  is  the 
reflection  coefficient  of  a  smooth  surface.  The  amount  of  power  in  the  diffused  field  may 
also  be  a  function  of  polarization.  The  values  of  Rd  that  were  measured  in  experiments 
is  within  the  range  of  0.2-0.4,  while  the  mean  values  were  between  0.3-0.35.  When 
6  =  0  (grazing  incidence),  A</»  =  0,  according  to  Equation  3.21;  and  <  R2S  >=  1. 
Under  these  conditions,  p8  =  Dp0.  This  means  that  the  specular  (coherent)  component 
is  equal  to  the  value  of  the  field  for  a  smooth  surface.  As  S  increases,  the  diffuse 
(incoherent)  component  increases,  until  the  phase  difference  approaches  the  value  of 
120  mradians.  From  this  point  the  diffused  component  decreases.  Also,  an  increase  in 
S  decreases  the  specular  component  (see  Equations  3.19-3.21).  [Ref.  2,4] 

For  sea  surface,  one  of  the  models  for  forward  scattering  expresses  the  electric 
scattered  field  as  a  vectorial  summation  of  the  specular  and  diffused  fields  with  zero 
mean.  The  phase  of  the  specular  field  is  equal  to  the  phase  for  a  smooth  sea;  its 
magnitude  is  reduced  by  the  scattering  coefficient  (Equation  3.19).  [Ref  3] 

Measurements  were  made  by  Beard  [Ref.  12]  on  the  reflected  field  above  sea 
surface,  for  roughness  varying  from  smooth  sea  — y*-  <  |  till  the  value  of  — jp-  = 
^.  The  transmitting  wavelengths  were  0.86,  4.2,  and  5.3  cm.  Both  polarizations, 
vertical  and  horizontal,  were  introduced  separately.  The  measured  results  indicated 
that  [Ref.  2]: 
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Figure  3.8:  Dependence  of  Mean  Square  of  Scattering  Coefficient  on  Phase 
Variance,  for  Specular  Reflection  from  Rough  Surface  [Ref.  2]. 
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1)  The  maximum  values  of  power  contained  in  the  diffused  field  is  approximately 
30%  of  the  specular  reflection. 

2)  For  both  polarizations,  the  results  did  not  change  the  polarization  independence. 
Probably  the  most  common  distribution  of  forward  scattered  field  from  a  sea 

surface  is  the  Ricean  distribution;  which  is  a  result  of  a  model  where  the  amplitude 
fluctuation  of  the  forward  scattered  diffused  field  is  Rayleigh  distributed,  and  its  phase 
is  uniformly  distributed.  This  model  has  proven  to  be  a  good  approximation,  also,  by 
experimental  results.  It  should  be  noted  that  for  very  low  grazing  angles,  the  ratio  of 
the  diffused  to  specular  powers  is  2-3  time  larger  than  for  higher  incidence  or  grazing 
angle.  [Ref.  3,6] 

The  root  mean  square  (RMS)  of  the  sea  wave  height  Ah  can  be  deduced  by 
observation  or  from  sea  state  tables.  These  values  are  given  in  Table  3.1,  and  will  be 
used  in  Chapter  IV  to  calculate  the  coverage  pattern  of  the  radar  due  to  multipath. 

TABLE  3.1:  PARAMETERS  OF  SEA  STATE.  [REF.  1,7]. 


Sea 

Type  of 

Wave 

RMS  Height 

State 

Roughness 

Height  [m] 

[cm] 

0 

calm 

0 

0 

1 

smooth 

0-0.3 

0-6.5 

2 

slight 

0.3-1 

6.5-21 

3 

moderate 

1-1.5 

21-32 

4 

rough 

1.5-2.5 

32-54 

5 

very  rough 

2.4-4 

54-86 

6 

high 

4-6 

0.86-1.3 

G.     SHADOWING  EFFECT 

When  radar  operates  at  very  low  grazing  angle  (<  1°),  the  scattering  from  the 
surface  will  be  "partial" .  This  means  that  large  parts  of  the  electromagnetic  wave  lie  in 
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the  geometric  shadow  of  the  sea  wave  crest  [Ref .  3] .  The  trough  of  the  wave  is  not  seen 
directly,  though  some  of  the  shadowed  areas  are  seen  due  to  refraction  illumination. 

However,  for  grazing  angles  less  than  1°  the  shadow  effect  greatly  affects  the 
backscatter  echo. 

Two  analytical  approaches  were  developed  for  the  shadowing  problem.  A  statisti- 
cal exercise  in  geometrical  optics  and  a  diffraction  problem  in  which  geometrical  optics 
shadowing  theory  is  considered.  The  procedure  of  the  statistical  approach  involves 
many  statistical  expressions,  but  eventually  results  in  a  shadowing  function  which  is 
the  percentage  of  the  illuminated  surface  area  at  a  specific  grazing  angle.  At  small 
angles  it  is  a  linear  function  of  the  angle.  This  approach  has  not  yet  been  proven  as 
representing  reliable  practical  shadowing  effects  at  low  grazing  angles.  [Ref.  6] 

In  the  geometrical  shadowing  theory  mentioned  above,  it  was  found  that  at  ex- 
tremely low  grazing  angles  (<  1°)  the  sea  surface  resembles  a  dark  plane  with  illumi- 
nated "spots".  The  size  and  location  of  these  spots,  within  the  reflecting  plane,  are  a 
function  of  the  wind  speed  and  the  grazing  angle.  This  result  is  closer  to  the  practical 
radar  picture,  involving  such  phenomenon  as  sea  spikes,  rather  than  the  shadowing 
function  picture  which  presents  diffuse,  vanishing,  scatter.  However,  this  geometrical 
theory  should  be  checked  carefully  before  applying  it  to  practical  cases.  [Ref.  2,6] 

H.     DEPOLARIZATION 

Depolarization  is  defined  as  a  change  in  the  polarization  of  an  electromagnetic 
wave  due  to  propagation  (in  a  medium  other  than  free  space)  and  including  surface 
diffraction,  reflection,  and  scattering. 

In  describing  quantitatively  the  depolarization  scattered  from  an  object,  a  so- 
called  scattering  matrix  is  derived.  The  scattering  matrix  describes  completely  the 
phase  and  amplitude  of  the  reflected  electromagnetic  wave  at  the  receiver.  For  better 
understanding,  suppose  that  an  electromagnetic  plane  wave  is  transmitted  with  its 
electric  field  in  the  y  direction  (see  Fig  3.9),  and  represented  by  E^.  Let  the  reflected 
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field  from  an  object  (or  a  surface)  at  the  receiver,  be  represented  by  E*,  and  be  given 
by  [Ref.  8]: 


Erx  =  PyxEl,  (3.25) 

where  pyx  represents  the  reflection  coefficient  of  the  incident  surface  and/or  object. 
(The  first  subscript  stands  for  the  transmitted  field  and  the  second  subscript  stands  for 
the  receiving  component.)  The  reflection  coefficient  pyx  is  complex,  and  includes  the 
change  in  phase  and  amplitude,  of  the  incident  field. 

Since  any  polarization  of  a  plane  wave  can  be  described  by  two  orthogonal  com- 
ponents, polarized  linearly,  a  scattering  process  can  be  described  completely  by  the 
polarization  scattering  matrix  [after  Ref.  8]: 


v  J 


(3.26) 


*^jb  Pxx     Pvx 

El\  I  Pry      Pyy 

Figure  3.9  illustrates  two  linear  and  orthogonal  electric  fields,  which  produce  by  sum- 
mation a  right  circular  electric  field.  For  a  perfect  conducting  plane  surface,  pxx  should 
be  equal  to  pyy.  Also,  the  scattered  energy  (square  electric  field)  from  a  surface  is 
mainly  a  function  of  the  transmitted  wavelength  and  the  reflecting  surface  orientation 
and  roughness.  [Ref.  8] 

In  1963  it  was  shown  by  Beckman  and  Spizzichino  [Ref.  1]  that,  if  a  plane  wave  is 
purely  polarized  (horizontally  or  vertically),  no  depolarization  will  occur  in  the  reflected 
wave  in  the  plane  of  incidence.  If  the  same  incident  wave  will  be  scattered  out  of  the 
incident  plane,  it  will  be  strongly  depolarized.  These  scattered  components  are  a 
result  of  the  irregularities  of  the  incident  surface,  and  related  to  the  diffused  scattering. 
[Ref.  1] 
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Figure  3.9:  Right  Circular  Polarization  Produced  by  Horizontal  and  Verti- 
cal Polarizations  [Ref.  8]. 
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According  to  experiments  above  sea  surface,  with  near  purely  polarized  wave,  the 
received  depolarized  specular  component  is  very  small  (at  least  25  dB  below  the  incident 
polarization).  On  the  other  hand,  the  diffused  scattered  field  is  strongly  depolarized. 

From  these  results  one  can  conclude  that  for  forward  scattering  from  sea  surface, 
depolarization,  for  linear  polarization  is  small.  This  fact  can  also  be  verified  when 
looking  at  echoes  reflected  from  targets,  especially  on  a  smooth  sea.  These  echoes 
might  fluctuate  sharply  due  to  interference  of  the  direct  and  the  reflected  wave.  The 
interference  can  occur  only  if  the  reflected  field  has  largely  the  same  polarization  and 
magnitude  as  to  the  direct  field.  It  is  shown  in  Ref.  1  that  very  little  depolarization 
takes  place  in  the  direction  of  specular  rays.  Therefore  the  interference  effect  is  addi- 
tional evidence  that  depolarization,  due  to  reflection  from  sea  surface  is  small  [Ref.  2]. 
This  phenomena  will  be  discussed  and  analyzed  in  Chapter  IV. 
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IV.  LOBING  PHENOMENA  ANALYSIS 

A.  INTRODUCTION 

In  Chapters  II  and  III  we  have  described  the  characteristics  of  earth's  surface  and 
a  little  about  air  refraction  (for  more  detail  see  Chapter  V),  which  dramatically  affects 
electromagnetic  propagation. 

In  this  chapter  we  will  discuss  and  analyze  the  influence  of  these  characteristics 
on  a  Naval  search  radar  (shipborne,  shore  installation  and  low-flying  aircraft). 

We  can  assume  that  the  detection  range  of  a  very  low  altitude  target  will  be 
limited  due  to  multipath  propagation  and  earth's  curvature  [Ref.  11]. 

In  this  chapter  we  will  discuss  and  analyze  the  lobing  phenomena,  caused  by 
multipath,  based  on  the  parameters  described  in  Chapter  II  and  III.  We  will  consider 
specifically  the  sea  surface,  and  the  detection  conditions  of  a  Naval  search  radar,  which 
usually  views  the  sea  at  extremely  low  grazing  angles  (<  1°).  The  covering  pattern 
of  detection  for  these  radars  is  strongly  affected  by  the  multipath  of  the  propagating 
electromagnetic  wave.  The  multipath  causes  constructive  and  destructive  interference 
due  to  earth  reflection.  This  effect,  together  with  media  and  radar  parameters,  will  be 
computed  to  illustrate  predicted  detection  patterns  of  the  radar. 

We  will  complete  the  analysis  with  some  practical  conclusions  and  recommenda- 
tions for  some  control  on  the  pattern  coverage  for  a  Naval  radar. 

B.  MULTIPATH  FADING  PHENOMENA 

The  main  sources  of  multipath  fading  are  reflections  from  the  earth's  surface  and 
from  refracting  layers  in  the  atmosphere.  In  the  commonly  known  multipath,  there  are 
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two  such  paths-a  direct  path  of  the  transmitting  wave  from  the  radar  to  the  target  and 
an  indirect  path  from  the  radar  to  the  target  via  reflection. 

The  geometry  of  this  scenario  is  shown  in  Figure  3.1.  From  this  figure  it  is  seen 
that  the  direct  and  indirect  wave  have  traveled  different  paths.  If  these  waves  have  the 
same  order  of  magnitude  they  can  strongly  interfere.  They  will  have  phase  differences 
caused  by  path  differences.  A  quantitative  analysis  of  the  lobing  phenomena  was  done 
by  Blake  and  is  described  herein  [after  Ref.  9]. 

Destructive  interference  will  occur  when  the  absolute  value  of  the  phase  difference 
is  between  |  +  27rn  and  y  +  27m  radians  (n  is  integer).  When  the  phase  difference  is 
7T  radians  (180°)  and  the  amplitudes  of  the  direct  and  indirect  field  are  almost  the  same, 
a  complete  cancellation  will  occur  at  the  receiving  point.  Constructive  interference  will 
occur  if  the  absolute  value  of  the  phase  difference  is  between  0  +  2nn  and  ir/2  +  27rn. 
When  there  is  no  phase  difference  {(f>  =  0),  the  constructive  interference  will  be  of 
maximum  value  due  to  arithmetic  summation  of  the  amplitudes. 

This  phenomena  can  be  seen,  if  we  look  at  the  field  at  a  fixed  height,  as  a  function 
of  the  range.  We  will  cross  through  maxima  and  minima  of  the  field  pattern  of  the 
radar,  caused  by  interference. 

The  goal  of  this  discussion  will  be  to  derive  procedures  and  equations  in  order  to 
predict  radar  range  performance  under  multipath  conditions.  We  will  derive  the  lobing 
factor  F;  considering  the  flatness  of  the  earth,  sea  roughness  and  polarization.  This 
was  discussed  in  Chapters  II  and  III.  For  simplification,  we  will  assume  that  the  target 
is  far  enough  so  that  the  usual  far-field  approximations  are  valid. 
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C.     THE  LOBING  FACTOR 

The  analytical  treatment  of  the  interference  phenomena  has  to  be  in  the  complex 
plane,  since  we  are  dealing  with  amplitudes  and  phases  of  the  waves.  Basically,  the 
phase  difference  of  the  direct  and  reflected  wave  is  dominated  by  three  factors: 

1)  the  geometric  path  length  difference  6  (see  Equations  2.6,  2.7), 

2)  the  phase  difference  of  the  reflected  wave  due  to  the  reflection  process,  and, 

3)  phase  difference  due  to  radiation  of  the  antenna  in  direct  and  indirect  directions. 
The  amplitude  difference  can  be  caused  by  the  following  reasons: 

1)  The  Naval  radar  with  narrow  beam  in  the  horizontal  plane  is  usually  not  uniform 
in  the  vertical  plane  which  can  cause  different  field  intensities  in  the  direct  and 
indirect  waves. 

2)  Both  waves  propagate  spherically,  resulting  in  an  inverse  square  law  power  reduc- 
tion. Since  the  path  of  the  waves  is  different  they  will  be  attenuated  or  spread 
differently. 

3)  The  reflected  wave  is  usually  attenuated  because  the  sea  surface  reflection  coef- 
ficient is  less  than  unity. 

1.  Applications 

The  spreading  loss  is  usually  neglected  for  both  waves  due  to  the  small  path 
difference  6  between  the  direct  and  reflected  wave.  The  major  effect  of  S  is  on  the  phase 
difference.  Most  of  the  Naval  surface  surveillance  radars  are  pulsed  radars.  The  effect 
of  interference  can  occur  only  if  6  <  ct,  where  r  is  the  pulse  width  and  c  is  the  speed 
of  light  in  free  space.  For  surface  and  low  flying  targets  and  moderate  height  of  radar 
antennas,  6  <C  ct  [Ref.  5,9].  This  will  result  in  susceptibility  to  the  multipath  effect. 
For  radars  which  use  pulse  compression  techniques  the  criterion  of  6  <  ct  is  for  the 
compressed  pulse. 

2.  Combined  Phase  Difference  a,  and  Detecting  Range  Rd 

The  magnitudes  of  the  reflection  coefficient  p,  and  the  pattern  factor  of  the 
antennas  /(7),  are  major  effects  on  the  amplitude  of  the  reflecting  wave.  The  elevation 
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angle  between  the  radar  antenna  and  the  target/receiver  is  7  (see  Figure  4.1).  The 
ratio  of  the  electric  intensity  in  the  direction  of  the  reflected  ray  to  the  intensity  of 
the  maximum  beam  direction  (boresight)  is  /(7r),  while  /fa)  is  the  same  ratio  in  the 
direction  of  the  direct  ray.  Both  |/(7)|  and  p  are  bounded  in  the  range  of  0  and  1. 

Now  we  can  derive  a  phasor  expression  for  the  direct  and  indirect  wave  fields. 
The  direct  wave  field  can  be  written  as: 

Ed  =  /fa)3>,  (4.1) 

where  Eo  is  the  phasor  value  of  the  free  space  field  in  the  beam  maximum  at  the  target 
range  Rj.  Similarly,  the  reflected  wave  field  is: 

Er  =  f(lT)E0-pe->a,  (4.2) 

where  a  is  the  total  phase  difference  caused  by  three  factors  mentioned  in  Section  C  of 
this  chapter,  and  p  is  the  magnitude  of  the  reflection  coefficient.  The  combined  field  at 
the  target  is: 

E  =  ET  +  Ed,  (4.3) 

Substituting  Equations  4.1  and  4.2  into  4.3  yields: 

E  =  E0[fM  +  fMpe-ja],  (4.4) 

the  lobing  factor  F  is  defined  as: 

F  =  \E/Eo\.  (4.5) 

Dividing  Equation   4.5  by  E0  yields: 

F  =  |/fa)  +  PfMe-ia\.  (4.6) 

If  we  assume  that  a  includes  the  phase  difference  fiT  —  fid  of  the  pattern  factor, 
then  by  factoring  out  fd  and  using  Euler's  theorem,  the  lobing  factor  F  is  written  as: 


F  =  fd  Vl  +  x2-\-2xcosa  ,  (4.7) 
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where 

*  =  ^,  (4-8) 

fr  and  fd  axe  scalar  magnitudes  of  the  direct  and  reflected  pattern  factors. 

We  can  represent  the  pattern  factor  f(jd)  and  f(jr)  in  the  form  of: 

/(7-)    =    /*-'*  (4.9) 

/(7r)    =    /re"iflr,  (4-10) 

where  /?<*  and  /?r  are  phase  angles  and  /<*  and  /r  are  the  magnitudes  of  the  pattern 
factors. 

As  mentioned,  a  is  the  total  phase  difference  between  the  direct  and  reflected 
waves  and  can  be  expressed  as: 

a  =  —  +  <f>  +  /3r-pd,  (4.11) 

where  the  first  term  in  a  is  caused  by  path  difference,  the  second  term  (<f>)  is  caused 
by  the  reflection  process,  and  the  third  term  (/?r  —  /3j)  is  the  phase  difference  of  the 
pattern  factors. 

The  detection  range  in  presence  of  multipath  and  atmospheric  loss  L  can  be 
calculated  by: 

*-§£•  (4-12) 

where  Ro  is  the  free  space  detection  range.  For  analyzing  the  lobing  effect  we  will 
propose  that  no  atmospheric  loss  occurs,  and  Equation  4.8  is  simplified  to: 

R  =  FRq.  (4.13) 

Equation  4.13  is  a  realistic  approximation  for  radar  frequencies,  since  the 
lobing  factor  is  dominant  when  multipath  occurs  and  affects  the  detection  range  more 
than  the  atmospheric  losses.  Since  0  <  F  <  2,  it  is  seen  that  the  range  Rd  under 
multipath  conditions  is  sometimes  greater  than  the  free  space  range,  and  sometimes 
smaller  than  the  free  space  range. 
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D.     PATH  DIFFERENCE  DERIVATION 

The  actual  antenna  pattern  located  at  A  in  Figure  4.1  can  be  found  by  treating 
it  as  a  two  source  interferometer  caused  by  the  real  antenna  and  its  image  from  the 
ground  at  C  [Ref.  5].  The  range  Rj  is  the  flat  earth  approximation  range.  If  we  assume 
that  Rd  ^>  h[  then  the  difference  between  the  direct  path  AE  and  the  indirect  path 
ADE  is: 

6  =  2/i'1sin0,  (4.14) 

Recall  from  Chapter  III,  Section  A,  h[  is  the  height  of  the  radar  antenna  and  6  is  the 

incident  (grazing)  and  reflected  angle  of  the  indirect  ray. 

For  small  grazing  angles,  which  are  typical  of  sea  surface  search  radars,  sin  9  can 

be  expressed  as: 

sin#  =  M  +  &,  (4.15) 

ixd 

where  h[  is  the  height  of  the  target.    By  assuming  that  h'2  ^>  h[  and  knowing  that 
Rd  ^>  h\  the  further  approximation  can  be  made: 

sinflw-^.  (4.16) 

Rd 

By  substituting  Equation  4.16  into  4.14,  the  path  difference  can  be  approximated  by: 

x  ~  2h'i  •  h2 


Rd 

2tt  . 


(4.17) 


By  multiplying  Equation  4.17  with  the  wave  number,  k  =  3r;  the  result  is  the  approx- 
imate but  practical  phase  difference  A<f>,  for  a  plane  reflecting  surface: 

A0«y— ^.  (4.18) 

The  elevation  angle  jd »  and  the  depression  angle  yri  can  be  calculated  by  geometric 
considerations: 

7,  =  sin"1  &ZM,  (4.19) 

Kd 


and 


*  = sin_1  irrw-  (4-2°) 


40 


©L---.$L-SL__.@ 


Figure  4.1:  Multipath  Geometry  over  a  Plane  Reflecting  Surface  [Ref.  9]. 
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Again  by  using  the  assumptions  h'2  >>  h[  and  R&  ^>  h\  we  conclude  that  : 

Id  =  7r  =  0.  (4.21) 

By  using  the  identity  of  Equation  4.21  in  Equation  4.14  there  results: 

6  =  2h[  sinfri).  (4.22) 

This  geometric  solution  of  the  multipath  effect  is  not  complete  since  we  have  to 
consider  also  the  reflection  coefficient  which  is  a  function  of  the  grazing  angle  8. 

E.     PRACTICAL  LOBE  FACTOR  DERIVATION 

The  solution  for  the  lobing  factor  is  achieved  by  substituting  Equation  4.17  into 
Equation  4.11;  this  result  for  a  is  substituted  into  Equation  4.7  together  with  the 
magnitudes  of  the  pattern  factors  and  the  reflection  coefficient.  [Ref.  3,9] 

A  special,  practical,  and  simplified  result  is  obtained  for  a  perfectly  smooth  sea 
and  horizontal  polarization.  For  this  case,  Figures  3.6  and  3.7  indicate  that  p  ~  1,  and 
<{>  ~  7r  rad,  independent  of  grazing  angle. 

For  surface  search  and  navigation  radars  the  main  beam  is  narrow  and  directed 
at  zero  elevation  angle,  and  the  vertical  beam  is  broad.  Under  these  conditions  fa  =  fa 
and  fi  —  fr  =  1  and 

a  =  —  +  7T.  (4.23) 

Substituting  this  simplification  into  Equation  4.7  yields: 

F=  |v/2  +  2cos(27r<5/A  +  7r)|.  (4.24) 

By  trigonometric  manipulation: 

F  =  2|sin(7r6/A)|.  (4.25) 

By  substituting  Equation  4.21  into  4.24  we  obtain,  finally,  an  especially  useful  result 
for  the  lobing  factor,  F: 

F  =  2\  sin  [{2Trh[  sin7(f)/A]  |.  (4.26) 
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F.     DISCUSSION 

Equation  4.26  allows  us  to  predict  the  relative  lobing  pattern  without  knowing 
the  range  and  altitude  of  the  target.  The  lobing  factor  F  depends  only  on  the  elevation 
angle  jd •  Figure  4.2  is  a  plot  of  the  lobing  factor  F  showing  the  vertical  plane  coverage 
affected  by  typical  multipath  from  a  smooth  sea.  The  antenna  height  is  10  meters  and 
the  transmitting  frequency  is  100  MHz. 

Note  that  the  maximum  range  Rd  is  twice  the  free  space  range  Rq  due  to  con- 
structive interference  since  the  maximum  value  of  F  is  2.  On  the  other  hand,  when  the 
target  is  located  at  other  elevation  angles  than  the  maximum  beam  value,  the  detec- 
tion range  can  be  less  than  the  maximum  free  space  range.  When  the  target  is  in  the 
direction  of  nulls,  no  detection  will  occur. 

As  seen  in  Figure  4.1,  variations  of  detection  range  relative  to  free  space  are 
bounded  between  the  two  extremes,  0  and  2,  and  are  a  function  of  the  elevation  angle. 
At  a  lobe  maximum,  we  have  [Ref.  5,9]: 

7mM(n)  =  sin-l(2n^)'A         n  =  0,1,2...      .  (4.27) 

The  minimum  occurs  when  the  same  argument  is  equal  to  multiples  of  27m.  The 
elevation  angles  at  which  minima  occurs  is  given  by: 

nX 
7min(n)  =  sin"1  —         n  =  0,1,2...      .  (4.28) 

Equations  4.27  and  4.28  indicate  that  the  maxima  and  minima  angles  depends 
only  on  the  wavelength  and  antenna  height.  This  is  true  only  under  the  following 
assumptions  [Ref.  9]: 

1)  The  pattern  factor  phase  difference  fid  —  fir  is  constant  at  all  elevation  angles. 

2)  The  phase  of  the  reflection  coefficient  is  also  constant. 

3)  The  reflecting  surface  is  plane  and  not  curved. 

These  assumptions  are  good  approximations  for  horizontal  polarization  and  an 
antenna  pattern  with  zero  elevation  beam.   If  we  take  into  account  also  the  reflection 
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Figure  4.2:  Vertical  Lobe  Pattern  Caused  by  a  Plane  Reflecting  Surface. 
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coefficient  and  antenna  pattern  factors,  it  affects  the  envelope  of  the  covering  factor, 
as  we  will  see  in  Section  G  of  this  chapter. 

The  free  space  radar  equation  for  the  received  signal  should  be  modified  by  the 
lobing  factor.  Due  to  free  space  loss  the  received  power  PT  of  a  radar  (two-way  loss)  is 
proportional  to: 

Pr  «  £.  (4.29) 

Doubling  the  free  space  detection  range  Ro  caused  by  constructive  interference  is  equiv- 
alent to  increase  the  transmitting  power  or  the  sensitivity  of  a  radar  16  times  or  equiv- 
alent^ by  12  dB. 

Equation  4.28  yields  that  the  first  minima  occurs  at  zero  elevation  angle.  This  is 
not  true  for  a  curved  surface  in  which  the  first  minimum  is  not  zero  and  its  elevation 
angle  is  not  zero.  Near  the  horizon  when  the  grazing  angle  approaches  zero,  the  lobing 
factor  is  very  low  causing  a  very  low  field  strength  at  this  region.  The  returned  signal 
is  dominated  by  the  spherical  nature  of  the  wavefront,  resulting  in  scattering  which 
can  be  taken  into  account  by  the  divergence  factor  D  [Ref.  3].  This  is  the  reason  why 
targets  at  very  low  grazing  angles  (low  altitude),  much  under  the  first  maximum  lobe 
angle  are  difficult  to  detect.  Typical  targets  with  low  grazing  angles  with  respect  to 
shipborne  radars  are  anti-ship  sea  skimming  missiles.  This  is  the  reason  for  the  effort 
being  made  trying  to  predict  and  control  the  location  and  shape  of  the  first  interfering 
lobe. 

The  flat  earth  equations  that  were  presented  in  this  chapter  are  good  approxima- 
tions for  relatively  low  antenna  height  such  as  shipborne  searching  radars  where  the 
reflecting  point  is  close  to  the  radar. 

G.     LOBING  ANALYSIS  FOR  TYPICAL  CASES 

We  want  to  now  analyze  plots  of  detection  patterns  affected  by  multipath  for 
radars  operating  above  sea  surface.  In  Figure  4.2  we  saw  a  detection  contour  of  a 
low  frequency  radar  (100  MHz)  with  horizontal  polarization  and  antenna  height  of  10 
meters  operating  above  a  perfectly  smooth  surface.  Figure  4.3  is  the  detection  contour 
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of  the  same  radar  except  that  the  height  of  the  antenna  is  doubled  to  20  meters.  This 
has  an  effect  of  lowering  the  angle  of  the  first  lobe  and  doubling  the  number  of  lobes 
nulls. 

Figure  4.4  is  a  detection  plot  of  a  radar  under  the  same  conditions  as  Figure  4.3 
except  the  frequency  is  now  10  times  higher  (1000  MHz).  This  effect  also  lowers  the  first 
beam  and  reduces  considerably  the  uncovered  space  between  the  lobes.  The  number 
of  lobes  and  nulls  are  increased  by  the  ratio  of  frequencies  (10  times).  Also,  according 
to  conservation  of  energy,  the  lobes  are  10  times  narrower.  One  can  conclude  that 
the  "penalty"  for  lowering  the  first  lobe  by  increasing  the  antenna  height  and/or  the 
frequency  is  that  the  number  of  nulls  is  increased  proportionally  (no  received  echo  in  null 
angles)  and  that  the  lobes  are  narrowed  proportionally.  For  shipborne  and  shorebased 
surface  search  radars,  increasing  the  frequency  and  the  antenna  height  is  recommended 
for  improved  detection  of  targets  with  low  grazing  angles  like  small  surface  combatants 
and  low  flying  aircraft  /missiles. 

If  the  polarization  of  the  transmitting  wave  is  vertical,  the  covering  pattern  will 
be  reduced.  This  is  particularly  true  for  the  lower  lobes,  due  to  reduced  values  of  the 
reflection  coefficient  p,  at  low  angles  (see  Figure  3.4).  Also,  the  phase  difference  is  not 
180°  as  it  was  in  the  case  of  horizontal  polarization  (see  Figure  3.5),  resulting  in  higher 
angles  of  maxima  and  minima.  [Ref.  9] 

In  Figure  4.5  we  include  also  the  effect  of  the  antenna  pattern  of  a  rectangular 
antenna  with  its  electric  field  uniformly  distributed.  The  pattern  is  derived  from  the 
Fourier  transform  of  the  rectangular  uniform  distribution  resulting  in  a  ^^  shape.  By 
comparing  with  Figure  4.3,  we  see  that  this  has  the  effect  of  modulating  the  basic 
pattern.  Also  the  higher  order  lobes  have  a  reduced  coverage. 

We  now  investigate  the  influence  of  the  scattering  coefficient  Rs,  considered  in 
Chapter  III,  on  the  radar  range.  The  effect  of  this  parameter  is  to  reduce  the  detection 
range  progressively  at  higher  grazing  angles.  This  is  because  <  R2a  >  decreases  expo- 
nentially with  increasing  grazing/elevation  angle  9.    Recall  from  Equations  3.20  and 
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Figure  4.3:   Vertical  Lobe  Pattern  for  the  Same  Parameters  as  Figure  4.2 
Except  Antenna  Height  Increased  to  20  m. 
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Figure  4.4:   Vertical  Lobe  Pattern  for  the  Same  Parameters  as  Figure  4.2 
Except  the  Frequency  is  Ten  Times  Higher  (1000  MHz). 
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3.21  that  the  mean  square  variation  in  Rs  is  given  by: 

<  R]  >=  e-2(2'rA'l8in*/A)\  (4.30) 

Figure  4.7  is  a  plot  with  same  parameters  as  in  Figure  4.6  but  for  a  higher 
roughness  of  sea  surface.  By  increasing  the  roughness  of  the  sea-mean  wave  height 
to  0.6  meters  instead  of  0.25  meters,  the  entire  pattern  is  considerably  reduced  due  to 
an  increase  of  the  diffuse  scattering  over  specular  coefficient. 

In  Figure  4.8  we  want  to  show  the  detection  contour  of  a  radar  with  appreciably 
higher  frequency  (1000  MHz  instead  of  100  MHz).  The  elevation  angle  is  also  reduced 
since  the  beam  width  is  proportional  to  ^,  where  D  is  the  physical  dimension  of  the 
antenna.  The  plot  is  carried  out  for  a  perfectly  smooth  sea  and  horizontal  polarization. 
We  can  see  the  effect  of  the  symmetrical  antenna  pattern  on  the  detection  contour. 
The  effects  of  higher  frequency  are  [Ref.  5,9]: 

1)  the  lobes  are  finer, 

2)  more  closely  spaced,  and, 

3)  the  first  lobe  is  at  lower  elevation  angle. 

The  effect  of  sea  roughness  at  high  frequencies  is  slightly  different  from  that  at 
low  frequencies.  Figure  4.9  uses  the  same  radar  as  in  Figure  4.8  except  mean  sea  wave 
height  is  changed  from  0  to  0.6  meters.  At  high  frequency  the  roughness  of  sea  has 
a  stronger  effect  especially  on  the  higher  angle  lobes.  An  additional  effect  which  is 
not  included  in  the  program  of  these  plots  is  that  nulls  are  partially  filled  in  this  case, 
especially  in  the  side  lobes.  The  multipath  effect  disappears  due  to  a  considerable 
decrease  in  the  specular  reflection  coefficient  ps  at  high  angles  [Ref.  10]. 

If  the  signal  is  spread  over  a  band  of  frequencies,  the  nulls  may  occur  around 
one  frequency  and  not  in  another  frequency,  resulting  in  a  decrease  in  fading  time.  By 
"spreading"  the  multipath  effect  over  a  wide  band  (whitening  it,  in  a  sense)  we  can 
reduce  the  effect  of  no  detection  at  nulls  which  occur  for  radars  operating  at  a  single 
frequency. 
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Figure  4.5:  Detection  Pattern  of  Rectangular  Antenna  with  Uniform  Field 
Distribution. 
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Figure  4.6:    Detection  Pattern  Same  Parameters  as  Figure  4.5  except  In- 
cluding the  Scattering  Coefficient  Effect. 
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Figure  4.7:  Detection  Pattern  Showing  the  Effect  of  Sea  Roughness. 
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Also  at  1000  MHz  if  we  change  the  polarizations  to  vertical,  the  lobes  magnitudes 
would  reduce  at  elevation  angles  higher  than  2°  as  can  be  inferred  from  the  pv  curve 
in  Figure  3.4. 

Some  typical  Naval  search  radars  operate  in  the  S  band  (2-4  GHz).  Figure  4.10 
illustrates  this  practical  detection  pattern  at  ideal  conditions  for  lobing  phenomena 
over  a  perfectly  smooth  sea. 

A  very  useful  way  to  analyze  the  multipath  effects  on  the  detection  capability  of 
a  radar  is  a  plot  of  the  received  signal  versus  range.  This  method  is  especially  used  to 
analyze  targets  flying  at  constant  altitude  towards  a  radar.  This  situation  can  simulate 
the  detectability  of  an  anti-ship  missile  approaching  a  ship. 

This  plot  is  more  convenient  for  analysis  of  detection  capability  of  targets  in  the 
first  lobe  (very  low  altitude  targets)  within  the  horizon  range.  In  this  plot  —93  dBm 
is  the  minimum  detectable  signal  (MDS).  Signals  below  this  level  will  not  be  detected. 
We  see  that  the  "dips"  in  echo  signal  due  to  multipath  are  narrower  and  more  frequent 
as  range  decreases.  This  effect  has  been  experimentally  verified. 

We  will  only  mention  here  that  surface  reflection  causes  errors  in  height  measure- 
ments. This  is  because  the  reflected  wave  can  be  considered  by  the  radar  as  an  echo 
signal  from  a  target  located  at  the  same  range  as  the  actual  target  but  at  the  opposite 
height,  i.e.,  the  image  of  the  actual  target.  When  constructive  interference  occurs,  the 
radar  will  measure  the  mean  value  of  the  elevation  angles  of  both  the  actual  and  the 
image  targets,  resulting  in  a  severe  error  in  height  measurement. 
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Figure   4.8:     Detection   Pattern   for   Higher   Frequency   (1000   MHz)    and 
Narrower  Vertical  Beam  Width. 
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Figure  4.9:  Effect  of  Sea  Roughness  on  a  1000  MHz  Detection  Pattern  of  a 
Radar. 
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Figure   4.10:     Detection   Pattern   under  Lobing  Effect    for  S    Band   Naval 
Surface  Searching  Radar. 
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Figure  4.11:    Received  Signal  Relative  to  Minimum  Detectable  Signal,  for 
Multipath  and  Free  Space  Conditions. 
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H.     SUMMARY  AND  CONSEQUENCES 

We  have  seen  in  this  chapter  that  the  specularly  reflected  field  may,  if  it  is  com- 
parable in  amplitude  to  the  direct  propagating  field,  produce  deep  fades  in  the  total 
field  through  interference.  The  diffused  field  has  much  less  effect  on  the  fading  when  it 
occurs. 

Atmospheric  refractive  layers  can  cause  many  different  path  rays.  In  this  chapter 
we  considered  only  the  interference  caused  by  the  direct  ray  reflected  from  sea  surface, 
with  no  effect  of  the  atmospheric  refractive  layer. 

The  main  parameters  affecting  the  detection  pattern  of  a  radar  are: 

1)  Antenna  height:  Increasing  the  antenna  height  increases  the  number  of  lobes 
with  less  uncovered  spaces  in  the  detection  pattern.  The  first  beam  is  lowered 
resulting  in  better  detection  of  low  altitude  targets. 

2)  Transmitting  frequency:  Higher  propagating  frequency  will  produce  a  covering 
pattern  with  finer  and  more  lobes.  The  lowest  lobe  will  be  considerably  lower. 

3)  Polarization:  The  lobing  phenomenon  is  more  prominent  for  the  horizontal  po- 
larization since  the  magnitude  and  phase  of  the  reflected  waves  are  more  closely 
correlated  to  those  of  the  incident  wave.  For  vertical  polarization  the  reflection 
coefficient  (magnitude  and  phase)  depends  strongly  on  the  grazing  angle  and 
decreases  at  higher  angles  and  frequencies,  resulting  in  a  decreased  detection 
pattern.  The  effect  is  stronger  at  higher  frequencies  (above  1000  MHz). 

4)  How  to  decrease  the  lobing  phenomena?  For  some  operational  applications,  it  is 
possible  to  overcome  the  multipath  effect.  For  example,  in  the  detection  of  low 
altitude  targets  and  in  active  jamming,  fading  is  not  desirable.  Some  means  of 
reducing  the  fading  phenomenon  are: 

a^e  of  spread  spectrum  techniques  (such  as  frequency  agility  or  hopping),  result- 
ing in  a  spread  of  fading  time  over  the  spectrum. 

fc^ansmit  vertical  polarization  at  frequencies  above  1000  MHz.     At  very  high 
frequencies  (35  GHz),  the  multipath  effect  is  minimal. 

The  above  parameters  can  be  easily  controlled  by  the  radar  designer. 
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V.  ANOMALOUS  ATMOSPHERIC 
PROPAGATION 

A.  INTRODUCTION 

In  the  previous  chapters,  we  considered  and  discussed  mainly  the  characteristics 
and  effects  of  earth's  surface,  i.e.,  sea  surface,  assuming  normal  atmospheric  conditions. 

In  this  chapter  we  will  consider  propagation  under  abnormal  conditions  called 
ducting  or  anomalous  propagation  caused  by  rapid  decrease  of  the  refraction  index, 
especially  in  warm  areas  such  as  the  Mediterranean  Sea.  This  phenomenon  considerably 
extends  the  detection  range  of  a  radar  within  the  ducting  layer  compared  to  horizontal 
free  space  range.  On  the  other  hand,  from  conservation  of  energy,  the  increase  range  in 
one  path  will  cause  a  decreased  range  in  another  path,  resulting  in  holes  of  detection 
outside  the  ducting  layer. 

We  will  consider  mainly  the  surface  or  evaporation  duct  because  it  is  the  most 
likely  one  above  a  sea  surface.  Refraction  losses  are  also  discussed.  We  will  end  with 
some  conclusions  for  radars  operating  in  a  duct. 

B.  TROPOSPHERIC  REFRACTION 

Electromagnetic  waves  propagating  in  the  troposphere  do  not  travel  in  straight 
lines;  instead  they  tend  to  bend,  due  to  refractivity  of  the  troposphere,  which  is  related 
to  its  dielectric  constant.  This  dielectric  constant  is  a  function  of  temperature,  pressure, 
and  water  vapor  content  of  the  air.  [Ref.  9] 

The  two  main  effects  of  this  phenomenon  on  radar  performance  are: 

1)  An  increase  in  the  maximum  range  of  propagation  is  increased,  and, 

2)  errors  in  measurements  of  elevation  angle. 
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Figure  5.1  illustrates  a  simplified  scheme  of  these  effects. 
1.      Refraction  Effect 

Refraction  or  bending  of  electromagnetic  waves  in  the  atmosphere  is  caused 
by  the  variation  of  the  velocity  of  propagation  with  height. 

The  index  of  refraction  is  the  ratio  of  the  velocity  of  propagation  in  free 
space  to  that  in  an  actual  medium.  For  microwave  frequencies  propagating  in  the 
troposphere,  which  contains  water  vapor,  the  refractivity  can  be  calculated  by  [Ref.  5]: 

xr      /        ,n-„h      77.6p      3.73xl05e 

tf  =  (n  -  1)106  =  -^  + ,  (5.1) 

where  p  is  the  barometric  pressure  in  mbars,  e  is  the  partial  pressure  for  water  vapor 
in  mbars  and  T  is  the  absolute  temperature  in  degrees  K.  The  refractivity  N  is  used 
instead  of  the  index  of  refraction  n  because  it  is  a  more  convenient  unit  when  dealing 
with  propagation.  The  index  of  refraction  normally  decreases  with  height  since  the 
barometric  pressure  p  and  the  water  vapor  content  e  decrease  much  faster  with  height 
than  the  temperature  T. 

The  index  of  refraction  is  also  given  by: 

n  =  -,  (5.2) 

V9 

where  c  is  the  velocity  of  light  in  a  vacuum  (3  x  108  m/sec)  and  vg  is  the  velocity 
of  the  propagating  wave  in  the  actual  medium.  Near  the  surface  of  earth  the  index 
of  refraction  is  1.0003.  Since  n  decreases  with  altitude  the  velocity  of  propagation 
increases  with  altitude.  According  to  Snell's  law  the  propagating  wave  tends  to  bend 
towards  the  lower  velocity.  The  result  is  an  increase  in  propagating  range. 

Since  the  refraction  index  n,  varies  with  height,  it  can  be  shown  from  Snell's 
law  for  spherical  geometry  that  the  factor  k  is  varying  with  height  and  has  the  form  of 
[Ref.  3,5]: 


*  =  1-l     l A    7m\>  (5-3) 

1  +  re{dn/dti) 
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Figure  5.1:   (a)  Extension  of  Radar  Horizon;  (b)  Elevation  Angular  Error 
[Ref.  5]. 
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Figure  5.2:  (a)  Bent  Beam  due  to  Refraction;  (b)  Straight  Line  Propagation 
with  Effective  Earth  Radius  4/3re  [Ref.  5]. 
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where  dn/dh  is  the  gradient  of  the  index  of  refraction  with  height. 

In  Chapter  III  we  considered  the  refraction  effect  by  increasing  the  earth's 
radius  by  a  constant  factor  of  4/3,  to  obtain  the  effective  earth  radius  for  straight  line 
propagation.  This  is  only  an  average  radius  and  should  not  be  used  for  anything  other 
than  general  calculations.  It  gives  a  rough  idea  of  refractivity  under  normal  atmospheric 
conditions.  Figure  5.2  shows  the  replacement  of  earth  radius  by  an  effective  radius  of 
4/3re  assuming  a  homogeneous  atmosphere  and  constant  gradient  of  index  of  refraction 
with  height. 

The  variation  of  index  of  refraction  with  height  was  found  to  be  an  exponen- 
tial function  rather  than  a  linear  function  as  assumed  for  k  =  4/3.  The  exponential 
model  for  refractivity  is  given  by  [after  Ref.  5]: 

Na  =  Ne  exp  [-Ce(hT  -  hR)] ,  (5.4) 

where  Ne  is  the  refractivity  at  the  surface  of  the  earth,  Ce  is  a  constant,  Hr  is  the  radar 
antenna  height  and  hj  is  the  target  altitude.  Use  of  this  refractivity  model  results  in 
smaller  errors  in  elevation  angle  measurements  of  the  target. 

C.     REFRACTION  LOSS 

At  normal  refraction  conditions,  when  the  vertical  gradient  of  refractivity  is  con- 
stant (jjr  =  const.),  the  main  effect  will  be  an  increased  elevation  error  as  measured 
by  the  radar.  However,  another  effect,  which  is  considered  by  Weil  1973  [Ref.  9],  is  a 
divergence  loss  of  signal  strength  that  is  caused  by  the  refractions  of  rays  from  atmo- 
spheric inhomogenities.  This  loss  is  also  known  as  Lens  Effect  Loss  since  the  effects 
caused  by  the  atmosphere  are  similar  to  the  lens  effect. 

To  illustrate  this  effect  consider  a  source  radiating  four  rays  A,B,C,D  as  in 
Figure  5.3.  Without  any  absorption  or  refraction  disturbancies  at  the  distance  R,  these 
four  rays  will  be  spread  out  through  an  area  equal  to  6  R2,  where  6  is  the  angular 
separation  in  radians  between  the  rays.  Now  assume  that  the  rays  are  propagating 
under  refractive  conditions  which  will  cause  them  to  bend  rather  than  to  propagate 
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in  straight  lines.  The  rays  C  and  D  have  a  slightly  smaller  angle  with  respect  to  the 
vertical  gradient  of  refraction,  thereby  causing  them  to  bend  slightly  more  than  rays  A 
and  B.  This  will  result  in  a  larger  spreading  area  of  these  rays,  causing  a  decrease  in 
power  density.  [Ref.  9] 


WAVEFRON 


RAYS 


SOURCE 


Figure  5.3:  Ray  Pattern  for  Refractive  Loss  Analysis  [Ref.  9]. 


This  refractive  or  lens  effect  loss,  can  be  calculated  by  the  ratio  of  wave  front 


areas: 


A' 
L'  =  — 
L        A' 


(5.5) 


where  A'  is  the  wavefront  area  under  refractive  conditions  and  ^4  is  the  free  space 
wavefront  area. 

Figure  5.4  shows  the  lens  effect  loss  in  dB  as  a  function  of  the  range  of  propagation 
with  the  elevation  angle  as  a  parameter.  The  conclusion  from  this  graph  is  that, 
although  the  lens  effect  is  not  large,  it  is  important  in  the  case  of  long  range  detection 
radars  trying  to  detect  targets  at  low  altitudes  (low  angles).  At  angles  higher  than  2°, 
the  lens  effect  loss  is  seen  to  be  negligible. 
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Figure  5.4:    Lens  Effect  Loss  as  a  Function  of  Elevation  Angle  vs.    Range 
[Ref.  4]. 

D.     PROPAGATION  UNDER  ANOMALOUS  CONDITIONS 

Theoretically,  under  ideal  conditions  with  no  attenuation  and  absorption,  the 
vertical  gradient  of  index  of  refraction  could  have  a  high  value  causing  rays  to  curve 
similar  to  the  earth  curvature  resulting  in  continued  propagation  around  the  earth. 
The  effective  radius  will  be  almost  infinite  if  jjjr  =  —1.57  •  10"' m-1  allowing  horizontal 
rays  to  propagate  around  the  earths  curvature  with  no  horizon  limits.  This  condition 
is  called  "ducting,"  in  which  the  electromagnetic  ray  is  "trapped."  These  conditions 
can  occur  only  within  an  atmospheric  layer  with  a  finite  thickness. 

Values  of  vertical  refractive  gradient  for  normal  and  abnormal  propagation  are 
known  to  be  [Ref.  5]: 
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1)   For  normal  atmospheric  refraction  0  <  jr  <  —0.787  -10  7m  *  or  1  <  K  <  2 


dh 

2)   For  superefraction  conditions  -0.787  •  10"7  <  ^  <  -1.57  •  lO^m"1  or 
2<  K  <  oo. 

Ducting  occurs  when  4?  >  —1.57-  10-7m_1.  This  phenomenon  can  be  compared 
to  a  leaky  waveguide.  There  are  places  along  the  layer  in  which  the  index  of  refraction 
increases  with  altitude  instead  of  decreasing,  causing  the  ray  to  bend  upwards  instead 
of  downwards.  In  spite  of  the  "leakage,"  ducting  usually  increases  tremendously  the 
horizontal  detection  range  of  radars. 

Once  duct  conditions  are  fulfilled,  the  refractive  index  within  the  duct  is  greater 
than  the  threshold  for  ducting,  resulting  in  the  ray  curving  down  in  a  curvature  larger 
than  that  of  earth's  surface.  If  the  surface  is  water,  some  specular  reflection  will  occur. 
Since  the  ray  is  trapped,  this  process  will  proceed  continuously  within  the  duct  as 
shown  in  Figure  5.5.  It  is  shown  in  the  figure  that  there  is  a  critical  elevation  angle 
at  which  the  ray  propagates  within  the  ducting  layer.  Beyond  this  elevation  angle,  the 
rays  will  propagate  out  of  the  layer.  This  critical  angle  is  usually  less  than  0.5°  meaning 
that  the  rays  are  almost  parallel  to  the  surface.  [Ref  9] 

This  type  of  duct  with  a  low  layer  above  surface  is  called  a  surface  duct,  or 
evaporation  duct,  if  it  occurs  above  water  surface.  A  duct  which  occurs  at  higher 
layers  is  called  an  elevated  duct.  [Ref.  5,9,13] 

The  evaporation  duct  usually  occurs  within  a  few  tens  to  almost  a  hundred  feet 
above  the  sea  surface.  The  duct  effect  can  be  interpreted  as  "forcing"  the  wave  to 
propagate  in  one  dimension  (horizontal)  rather  than  in  the  two  dimensions  under  nor- 
mal conditions  (horizontal  and  vertical).  This  results  in  power  density  spreading  as  the 
inverse  first  power  of  the  range  (l/R)  instead  of  1/R2.  In  two  way  propagation  radar, 
the  receiving  energy  will  be  proportional  to  l/R2  under  ducting  conditions  instead  of 
the  l/R4  in  free  space.  The  result  is  that  for  a  radar  and  a  target  within  the  ducting 
layer,  the  detection  range  will  be  much  beyond  the  "horizon"  range  of  free  space  (see 
Fig.  5,6).  [Ref.  9,  13] 
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Figure  5.5:  Ray  Propagation  Pattern  under  Ducting  Conditions  [Ref.  9], 
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Though  there  are  losses  caused  by  absorption  and  leakage,  the  ducting  phenomena 
has  an  effect  of  increasing  the  detection  range  within  the  layer.  This  extended  detection 
range  occurs  because  more  energy  is  directed  into  the  ducting  layer  by  refraction,  while 
less  energy  will  propagate  in  other  directions  causing  "holes"  or  poor  detection  in  areas 
outside  of  the  ducting  layer.  We  can  conclude  that  a  surface  or  evaporation  duct  favors 
detection  of  surface  or  sea  skimming  targets,  and  does  not  favor  the  detection  of  targets 
with  elevation  angles  greater  than  1°. 

E.     SURFACE  DUCT  CHARACTERISTICS  AND  ATMOSPHERIC  RE- 
LATIONSHIP 

As  mentioned  in  Chapter  V,  Section  A.l,  the  refraction  index  decreases  with 
height  in  normal  atmospheric  conditions.  Under  certain  conditions,  especially  in  sub- 
tropical regions,  there  is  a  warm  layer  over  a  colder  layer  which  means  that  temper- 
ature increases  with  height  instead  of  decreasing  with  height.  This,  together  with 
normal  rapid  decrease  of  pressure  and  evaporation,  results  in  a  more  rapid  decrease 
of  refraction  index  than  under  normal  conditions.  Such  ducts  are  common  in  tropical 
and  semitropical  areas  over  cold  sea  water  with  a  maximum  height  of  at  most,  a  few 
hundred  feet.  [Ref.  3,9] 

Another  type  of  duct,  the  evaporation  duct,  is  due  to  evaporation  of  moisture 
from  the  sea  surface.  This  duct  is  usually  lower,  20-60  ft  on  the  average,  and  is  more 
reliable  than  the  previously  discussed  duct.  This  type  of  duct  is  common  in  the  warm 
areas  of  the  globe  such  as  the  Mediterranean  Sea,  for  instance. 

The  maximum  wavelength  of  an  electromagnetic  wave  able  to  propagate  within 
a  duct  depends  on  its  thickness  d.  A  duct  could  be  thought  of  as  a  high  pass  filter 
allowing  only  the  higher  frequencies  (shorter  wavelength)  to  propagate  abnormally. 
The  maximum  wavelength  that  can  propagate  in  evaporated  duct  with  thickness  d  is 
given  by  [Ref  5] : 

Am„  =  2.5(-J),/2<P'2.  (5.6) 
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Figure  5.6:    Relative  Received  Signal  Under  Duct  and  Free  Space  Condi- 
tions. 
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For  an  X  band  radar  (A  =  3  cm)  the  required  thickness  for  ducting  to  occur,  would 
be  10  m.  Since  ducting  layers  are  not  very  thick,  ducts  will  usually  occur  for  higher 
microwave  frequencies. 

Compared  with  waveguides,  the  "cut-off  wavelength"  for  a  duct  is  not  so  strict. 
Modes  with  "harmonic  wavelength"  can  also  be  propagated  by  ducts.  [Ref.  5] 

F.     CONCLUSIONS  OF  ANOMALOUS  PROPAGATION 

Finally  the  following  points  should  be  considered  when  designing  and/or  operating 
a  radar  in  anomalous  propagating  conditions: 

1)  Evaporation  duct,  when  used  with  the  appropriate  height  of  antenna  and  wave- 
length, can  extend  the  detection  for  surface  or  low  flying  targets  considerably 
beyond  the  expected  range  with  normal  atmospheric  conditions. 

2)  The  extended  horizontal  detection  range  will  cause  reduced  range  or  holes  in 
detection  at  other  elevated  angles.  This  follows  from  conservation  of  energy. 
This  can  cause  severe  coverage  problems  for  air  surveillance  and  early  warning 
radars. 

3)  Usually  radar  is  designed  under  the  assumption  that  clutter  is  likely  only  in  the 
range  close  to  the  radar.  In  ducts,  clutter  is  likely  also  at  longer  ranges,  resulting 
in  a  signal-to-clutter  problem  for  radars  without  techniques  to  cancel  it  (such  as 
pulse  staggering  against  second  time  around  clutter.) 
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VI.  SUMMARY  AND  CONCLUSIONS 

A.  SUMMARY 

In  this  study  we  have  tried  to  quantify  the  propagation  characteristics  of  an 
electromagnetic  wave  through  material  and  media  boundaries.  Random  irregularities 
of  the  sea  surface,  the  so-called  sea  roughness,  have  been  investigated  and  included  in 
the  specular  and  diffused  reflection  coefficients.  The  curvature  of  the  earth's  surface 
is  taken  into  account  by  way  of  an  effective  radius.  These  parameters  were  used  to 
compute  the  detection  performance  contours  of  a  typical  microwave  search  radar.  The 
multipath  effect  appears  as  a  lobing  factor  on  the  detection  range.  Further,  the  effects  of 
anomalous  propagation  through  the  atmosphere  are  also  studied.  The  work  illustrates 
that  the  propagating  medium  has  a  significant  effect  on  the  detection  range. 

B.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  study  show  that  the  reflected  field  from  the  sea  surface  interferes 
with  the  direct  field,  and  depends  on  wavelength,  polarization,  sea  state,  angles  of 
incidence,  and  the  electric  properties  of  the  sea  surface. 

By  appropriate  design  and  operation  of  the  radar,  we  can  have  some  control  on 
the  detection  pattern. 

Frequency  and/or  spatial  diversity,  of  the  transmitting  frequency  and  antenna 
heights  respectively,  were  proved  to  be  very  effective  in  "filling  up"  dips  in  the  receiving 
signal  caused  by  multipath.  Also,  increasing  the  frequency  and  antenna  height  lowers 
the  first  lobe  of  the  detection  contour,  for  better  detection  of  targets  at  low  grazing 
angles. 


71 


In  areas  such  as  the  Mediterranean  and  the  Red  Seas,  abnormal  atmospheric 
conditions  for  propagation  are  very  common  for  most  parts  of  the  year.  In  order  to 
gain  the  benefit  of  the  ducting  phenomena  one  must  predict  the  thickness  of  the  ducting 
layer. 

By  using  appropriate  antenna  height  and  wavelength,  detection  of  targets  within 
the  ducting  layer  are  maximized.  Further  studies  should  be  made  for  predicting  ducts 
more  accurately  so  that  one  can  rely  on  this  phenomena. 
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APPENDIX:  DERIVATION  OF  ELLIPSES  FOR 

FRESNEL  ZONES 

To  derive  the  equations  of  the  ellipsoid  it  is  easier  to  change  coordinates  to  x\  y\  z' 
by  tilting  the  x  coordinate  by  the  angle  of  0  in  order  to  align  the  coordinate  with  R 
(see  Figure  2.3).  y'  will  be  parallel  to  y  and  z'  is  obviously  orthogonal  to  the  x'y'  plane 
[Ref.  3,4]. 

The  expressions  for  the  new  coordinate  systems  x'  y\z'  are  [Ref.  4]: 


x'    =     Ix-  -jcosfl+  Iz-  jsinfl,  (A.l) 

y'    =    y,  (A.2) 

z'    =    -(x-7-)sm6+  (z-  Zl*Z2)cos6.  (A.3) 

For  given  S  the  equation  of  the  ellipsoid  is: 

a2(y'2  +  z,2)  +  b2x,2  =  a2b2.  (A.4) 

where 

r  sec  6  +  6 

a    = ,  (A.5) 

b    =    ^VS2  +  2rSsecO.  (A.6) 


and 


tan0=^ — -.  (A.7) 

r 


We  are  interested  in  the  Fresnel  zones  in  the  xy  plane  transferring  the  ellipsoid  Equa- 
tion A.4  back  into  x,y,z  coordinates  with  z  =  0  in  the  transformation  equations,  we 
get  a  family  of  ellipsoids: 

a262(62cos20  +  a2sin20-c2>)        /        i  x2 

\2       ,.        2   .   2, L  =lx--r-p)    (b2  cos2  0  +  a2  sin2  0  +  a2y2),      (A.8) 

bl  cos2  0  +  a1  sin  0  \         2  / 
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where 


c    = 


P    = 


Z\  +Z2 


c(a2  —  b2)  sin  6  cos  0 


(A.9) 
(A.10) 


b2  cos2  B  +  a2  sin2  0 

Equation  2.5  represents  an  ellipse  corresponding  to  a  given  8,  the  center  of  the 
ellipse  is  on  x  axis  and  is  given  by: 


x0  = 


1- 


(£/r  +  sec0)2-l 


y0  =  o. 


By  setting  x  =  xq  the  semiminor  axis  is: 


Vi  =  ±&\  A 


2/3' 


b2  cos2  0  +  a2  sin^  6 


or 


yi  =  ±s 


2\ 


r  /  r  V  V      r 


1- 


(*i  +  *2)2 


^  [f-f^l+C^-)2]    -1 

By  setting  y  =  0  the  semimajor  axis  is  derived: 

/Tf 

Xi  =  X0  ±  t/i 


V^os^+T^sin2^ 
Substituting  Equations  A. 5,  A. 6,  and  A. 7  into  Equation  A. 14  yields: 


#1  —  Xq  ±  yi 


N 


i  + 


f  +  yi  +  (^L): 


-i 


(A.ll) 


(A.12) 


(A.13) 


(A.14) 


(A.15) 


For  practical  computations  these  equations  can  be  simplified  by  assuming  that 
Z\ ,  z2  and  the  path  difference  of  8  are  much  smaller  than  r.  This  yields  to: 


S+M*¥* 


28      fz2±z1  , 

-1«-+    — L  )    <  1. 

r        \      r 


(A.16) 
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